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PROBABILITY OF SATELLITE COLLISION 
SUMMARY 


A technique for computing the probability that a collision will occur 
between a given satellite (target satellite) and any one of the total population 
of debris satellites is presented. This method is restricted to near-circular 
target satellite orbits. The following assumptions were made: 

1. The satellite sample consists only of objects trackable by North 
American Air Defense Command (NORAD) . 

2. The satellite population remains static as of a certain epoch. 

3. Effects of aerodynamic drag, solar wind, and deliberate propulsive 
maneuvers are neglected. 

4. The target satellite and debris satellite ascending nodes are random 
over 360 deg. 

5. The arguments of perigee of the debris satellites are random over 
360 deg. 

6. The position of lapping between the target satellite and each debris 
satellite is random over 360 deg. 

7. The radius of apogee, radius of perigee, and inclination of each 
debris satellite and the target satellite remain constant. 

In order to make the probability calculations, each debris satellite was analyzed 
separately, and the total probability of collision was determined from the 
individual probabilities of each debris satellite. 

A realistic sample problem (the proposed modular Space Station) is 
presented to provide further insight into the collision problem. A satellite 
sample of September 1970 containing 1805 objects was used. For a Space 
Station mission of 10 yr and an inclination of 55 deg, the collision probability 
varied from approximately 0. 01 to 0. 04 over the altitude range being considered 
for the Space Station. 

The results of a parameter study to determine what types of debris 
orbits (size and shape) contribute a large part to the total collision probability 
are given. Near-circular orbits which cross the target vehicle altitude con- 
tribute the most to the total collision probability. 



INTRODUCTION 


Man’s understanding of Iris environment has advanced significantly 
with the advent of artificial earth satellites as remote sensing devices. When 
these instruments are launched into earth orbit, there is usually some debris 
(payload shrouds, explosive bolt fragments, spent stages, etc. ) which 
accompanies them. The rate at which these objects (payloads and debris) 
have been placed into earth orbit has exceeded the rate at which they have 
fallen back to earth (due to orbital decay) , affecting a continual growth in the 
total number of artificial earth satellites. NORAD has identified 2588 objects 
in earth orbit as of October 31, 1971 | 1 | . If the present trend continues, 
this number will continue to grow. 

In view of this large and increasing satellite density in earth orbital 
space, the question arises as to the probability of a collision between two of 
the satellites. This question has particular significance for manned missions 
since crew safety is always a high priority item. The protection of expensive, 
unmanned satellites (communications, meterological, etc.) against collision 
with other objects would also be desirable. Decisions with regard to M safe M 
altitudes or necessity of collision avoidance systems might be intimately 
dependent on a realistic assessment of the collision hazard. 

This report presents a method of determining the probability that a 
given satellite (target satellite) will collide with any other satellite (debris 
satellite) in earth orbit. The computational technique is described in general 
terms in the text, and a realistic sample problem (the proposed modular 
Space Station) is analyzed to provide further insight into the collision problem. 
The equations are derived in the seven appendices. 


PROBLEM DESCRIPTION 


The method of computing collision probability for a single debris 
satellite with a target satellite is first shown, and is then extended to include 
the entire debris satellite population. 

Figure 1 shows an impending collision between the target satellite and 
a debris satellite. Clearly, a collision possibility exists only in the vicinity 
where the paths of the two bodies intersect (within A'k /2 as shown) . This 
implies that both bodies must be close to the line of intersection between the 
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orbit planes and that their altitudes be nearly the same. Further, both 
bodies must be there at nearly the same time. A collision is possible only 
at such times that one body is in the process of lapping the other body (lapping 
is the acquisition of a zero phase angle between the objects) . 



Figure 1. Orbital geometry for a collision 
situation between two orbiting objects. 
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In this study, the target satellite is simulated as a sphere of radius R. 
The dimensions of the debris satellite are neglected explicitly, but can be 
accounted for in an actual case by increasing the radius of the target sphere 
by an appropriate amount. If the debris satellite touches or intrudes into the 
sphere representing the target satellite, then a collision is defined. The 
immediate problem, then, is to determine the probability that this will occur. 


MATHEMATICAL FORMULATION 


For this problem, it is convenient to define a coordinate system at 
the center of the target satellite as shown in Figure 2. This coordinate 
system will be referred to as the ’’relative” coordinate system. The y axis 
of this system remains in the direction of the position vector R^, and the 

z axis is in the direction of the momentum vector R^ x V . The x axis is 

a curvilinear axis which is tangent to the orbit path and remains parallel to 
the local horizontal. Defining the relative coordinate system in this way 
restricts this study to circular orbits for the target vehicle. 



For specified orbits and a given mutual orbit plane orientation (specific 
values of inclinations, arguments of perigee, and ascending node separation) , 
the velocity of the debris satellite as it passes through the orbit plane of the 
sphere can be determined in the relative coordinate system (Appendix A) . 

If the path of the debris satellite is approximated as a straight line in the 
vicinity where collision is possible, then the direction of this relative velocity 
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vector represents the direction of the trajectory of the debris object* When 
this trajectory passes through the target sphere, a collision occurs. 


The debris satellite lias an opportunity to pass through the sphere only 

when lapping between the two objects is occurring. This is made clearer by 

examining Fig. 3 which shows an inertial view of the paths of the debris and 

target satellites., Assume that the target satellite is in position 1 when the 

debris satellite passes through the target orbit plane. IT the period of the 

target satellite is less than the period of the debris satellite, thru the target 

satellite will have traversed a complete orbit plus some distance AS (position 

2) at Ihe time the debris satellite passes through, the target plant' again. 

Therefore, with each successive orbit of the debris satellite, the target 

satellite approaches AS closer to the position in its orbit Lit which a collision 

can occur. This same idea, is illustrated in Fig. 1 with respect lo the relative 

coordinate system, This figure shows the debris satellite passing through the 

target orbit n'sne at an altitude difference of y _ . With each orbital period, 

1 ' * d 

the debris satellite passes through the target plane AS closer to the sphere. 


In this relative reference frame (Fig* 4) there 


ter is ties common 


each debris satellite pass through 


arc* certain eharac- 
ihr target orbit plane: 


1. The value of y , (the altitude difference between die two bodies) is 
“ d 

approximately the same each pass. The reason fur this is Him the altitude 
of the debris object during passage is determined by Ihe position of the inter- 
section line in the debris orbit, which in turn is determined by the argument 


ot perigee (the angle between the ascending node and 
Perturbations in argument of perigee due to Ihe effec 
m ss are too slight to cause significant changes in y 


the line of perigee) . 

Is of the earth 1 s oblale- 
froui one pass to another. 


2. The direction in which the debris satellite is traveling as it passes 


through the target piano remains approximately the s 
passage. The same reason given for item 1 above? at 
a given angle between the orbit planes, the direction 


ame from passage to 
so applies here since, 
of travel is a. function 


for 

of 


the satellite T s flight path angle and the flight path angle during passage h 


function of argument of perigee. 


The distance between successive passes, AH, is a constant because, 
as stated earlier, this study is restricted to circular orbits of the target 
satellite. The path speed of the target satellite in circular orbits is a constant; 
thus, the net change in its position during one orbital period of the dTwjs 
oujoet will remain constant. 


D 




Net position change of the target satellite during 
one orbital period of the debris satellite. 
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4. The path of the debris satellite is almost a straight line in the 
vicinity of the target pjane where collision is possible (at least a distance 
2R from the target plane) . 

Assuming that the debris satellite’ s R 4 and R remain constant, the 

A p 

value of y d depends upon the orientation of the major axis of the debris satel- 
lite’ s elliptical orbit with respect to the line of intersection; i.e. , y depends 

d 

upon the ’’true anomaly" of the line of intersection as measured in the debris 

satellite orbit plane (Fig. 5) . Most orientations, as represented by the 

dashed line in Figure 5, will cause y , to be of such a value that it would be 

d 

impossible for the trajectory of the debris satellite to pass through the target 

sphere. However, when the intersection line occurs within some interval of 

true anomaly, At], approximately the position of equal altitudes (the shaded 

area) , a non-zero collision probability exists and can be calculated in the 

following manner. As shown in Figure 6, there will be a certain interval (d) 

along a line parallel to the relative x axis, such that if the trajectory of the 

debris object intersects the target orbit plane within this interval, the trajectory 

will pass through the target sphere. The exact value of d is dependent upon the 

value of y and the three-dimensional direction that the debris satellite is 
d 

traveling. The derivation of d is given in Appendix B. 

Once this collision distance ( d) has been calculated, the probability 
of collision can be calculated in the following manner. In Figure 6, the 
successive debris passes through the target orbit plane are shown occurring 
at particular positions along the relative x axis. The exact positioning of 
these passes is actually unknown and is therefore considered to be completely 
random. That is, the debris satellite passes can be shifted any amount up to 
a distance of AS (after which any subsequent shift produces positions which 
are a repeat of a previous position) , and each of these positions is equally 
probable to occur. 

To gain insight into how the probability of collision can be calculated 
from this, consider the following example. Assume that the interval AS is 
divided into six equal segments of length Ax (Fig. 7) and that there are a 
large number of cases (say 1000) of a debris satellite passing through the 
target plane in any random position within AS. This quality of randomness 
dictates that the same number of passages are possible in each Ax interval, 
because passage through any segment of AS is equally probable with passage 
through any other segment of equal length. In Figure 7, Ax= AS/6. Let 
m be the number of cases which can occur in the interval Ax out of the total 
number of cases, n, which occur within AS. Then 
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Equation (l) is for a Ax equal to ( l/6 ) AS* For any size Ax 


m 


n 

AS/Ax 


and 


m Ax 

n AS 


( 2 ) 

(3) 


A fundamental definition of probability is as follows: If there are n mutually 
exclusive, exhaustive, and equally likely cases, and m of these are favorable 
to an event A. then the probability of A is m/n. Since equation (3) holds for 
any Ax, it holds when Ax~- d. By the above definition, the probability that 
the debris satellite will pass within the interval d is 

P = HI _iL 

lap n AS 



Since the debris satellite has a single chance to pass through the interval d 


each time the debris and target satellites lap each other, 
bility per lap. 



is the proba- 



Figure 7. Geometrical explanation of collision probability expression. 


11 


Equation (4) was developed for a given mutual orbiv plane orientation, 

(the assumption was made that inclinations, arguments of perigee, and node 

separation angles were known) . The quantity, d, and therefore P, , is a 

lap 

function of these quantities. The perturbation rates (because of the earth' s 
oblateness) of node separation and argument of perigee are great enough so 
that it is impractical to maintain up-to-date measurements on these quantities. 
Therefore, for a given time (at which it is desired to compute a collision 
probability) , it is assumed that node separation A S 2 and argurrient of perigee 
co are unknown and that all values of these angles between 0 deg and 360 deg 
are equally probable. Rewriting equation (4) as functions of these variables 
yields 


lap 


d (co , Afi ) 
AS 


(5) 


For each value of A S 2 , only a very restricted range of values of co will yield 

a non-zero value for d; hence, for most values of co , P =0. The average 

lap 

collision probability per lap over node separation and argument of perigee is 

2jr 2t r 


lap 


1 

4tt 7 


f f P lED <^ d(ASi) 
0 0 la P 


( 6 ) 


An attempt was made to perform the double integration in equation ( 6 ) 
analytically. Since this attempt was not successful because of the complexity 
of the probability expression, the expression was numerically integrated by 
the trapezoidal method on a digital computer. Numerically integrating equa- 
tion ( 6 ) over both variables required an excessive amount of computer time 
because these calculations must be made for each debris satellite in earth 
orbit ( about 2400 presently) . This problem was circumvented ( losing very 
little accuracy) by noting, after several numerically integrated runs, that 
the plot of probability as a function of argument of perigee, co , manifests 
itself as four humps which resemble half ellipses (Fig. 8 ) . The reason for 
the four humps is made clearer by referring to Figure 5. Varying the argu- 
ment of perigee, co , from 0 deg to 360 deg has the effect of varying 17 , the 
orientation of the line of intersection between the two orbit planes relative 
to the position of the debris satellite perigee. The sum of co and 
17 is the argument of latitude (measured in the debris orbit plane) which 
remains constant for a given AS 2 . The shaded areas in Figure 5 are the 
range of intersection line positions where a non-zero collision probability 
exists. Each of these orientation ranges occurs twice when 77 (or co ) is 
varied from 0 deg to 360 deg when 
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and 


1 . r] is between ri and rj 
q r 


2. rj is between 180 + r\ and 180 + n . 

q r 

If the integral of this function (Fig. 8) is approximated as the area of 
two ellipses, then only the probability as a function of node separation need 
be numerically integrated. The method of determining the dimensions of 
the humps for the area calculations is described in Appendix C. A mathe- 
matical justification for the use of this ellipse approximation is given in 
Appendix D. 

A few very special cases (debris orbits) produce a probability versus 
cd curve having a shape similar to the one in Figure 9 and must be numeri- 
cally integrated. The elliptically shaped humps become distorted so that the 
area under one of the humps can no longer be approximated as the area of a 
half ellipse. These types of curves are generated when either the apogee or 
perigee of a debris satellite orbit is very near (within a few diameters of 
the target sphere) to the altitude of the target satellite. The humps are 
actually connected as in Figure 9-b when the altitude of apogee or perigee 
falls within the radius of the target sphere. Cases like these for which the 
elliptical approximation does not hold are numerically integrated over argu- 
ment of perigee. Very little additional computer time is required since only 
a few (usually less than 10) cases like these are encountered for any given 
target satellite altitude. * 


The probability that the debris satellite (the jth satellite) will collide 

with the target satellite at least once within any duration mission (P ) can 

v m 

be calculated once the probability for each lap in the mission, [P (equation 

tap 

6)] is known by using the following equation: 


P 

m 

j 


= 1 - (1 


5 ,ap> L 


( 7 ) 


where 


L = the total number of laps that occur between the target satellite 
and the debris satellite. 

The probability that the target satellite will collide with any one of the debris 
satellites is 
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P TOT “ 1 


-n 

j=i 


(i - p ) 
m. 

J 


where 


k = the total number of debris satellites in earth orbit. 
Equations (7) and (8) are derived in Appendices E and F, respectively. 


SAMPLE PROBLEM 


The problem of collision with debris satellites takes on special 
significance for the proposed modular Space Station because the expected 
mission duration of 10 yr and the large Space Station dimensions provide 
much greater exposure. An analysis of the collision hazard of the Space 
Station is included to illustrate the of the computational technique and to 
provide additional insight into the collision problem. 

The launch and operation of the modular Space Station are still some 
years away. The satellite population in existence at the time the mission is 
in progress will obviously determine the collision probability. However, future 
satellite populations are difficult to predict since the satellite launch rate 
depends on economic, political, and military factors. Therefore, an arbi- 
trary launch date and a static satellite population were assumed. The satel- 
lite sample used in this analysis was the population and distribution in exis- 
tence as of September 1970. This satellite sample contains only objects 
which are tracked by NORAD. The sample contains 1S05 objects, but a 
collision is possible only with those objects which cross the Space Station’ s 
altitude. 


RESULTS 


The total collision probability was calculated for a range of Space 
Station altitudes from 200 km (108 n. mi. ) to 1400 km (756 n. mi. ) consider- 
ing the entire satellite population. The results are shown in Figure 10. In 
generating these results, the Space Station inclination was held constant at 55 
deg and the mission duration was assumed to be 10 yr. 
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The Space Station was represented as a sphere with a radius of 
50 m. Although the actual size and configuration of the Space Station had 
not been finalized at the time of this study, the assumption was made, based 
on current design proposals, that a miss distance of 50 m from the center 
of the Space Station would represent either a collision or a very near miss. 
Therefore, the results in Figure 10 can be interpreted to be the probability 
that a debris satellite will intrude within a 100-m-diam sphere at least once 
during a 10-yr mission. A more detailed plot of the same data is given in 
Figure 11. In this figure the number of debris objects which cross the vari- 
ous altitudes is plotted. The general trend of the collision probability is 
to increase as the number of satellites increases. Figure 12 is a plot of 
collision probability versus Space Station inclination. The altitude of the 
Space Station was held constant at 500 km (270 n. mi. ). 

The collision probabilities presented in Figures 10, 11, and 12 
must be regarded as minimum values, since there is considerable debris in 
earth orbit which is not trackable by NORAD due to the small size or the 
low reflectivity of the object. 

A parameter study was performed to determine the types of debris 
orbits (size and shape) which contribute a large paid to the total probability. 
Figure 13 shows the effect of varying perigee for constant apogee altitudes 
and Figure 14 is the result of varying apogee for constant perigee altitudes. 
The inclination of each debris orbit was held constant at 0 deg. The collision 
probabilities in these figures are not total probabilities but are individual 
satellite probabilities. These curves show that the largest contributions to 
the total collision probability are made by debris whose orbits are nearly 
circular and near the Space Station altitude. In other words, the closer the 
apogee or the perigee approaches the Space Station altitude, the higher the 
probability that a collision will occur (neglecting the effect of inclination). 

Figure 15 shows an expanded view of the peak of one of the curves 
of Figure 13. The collision probability reaches a peak and diminishes to 
zero as the satellite perigee passes through the altitude blind of the sphere. 
The same trend is also present in Figure 14 as apogee passes through the 
altitude band of the sphere. 


Encountering some near-circular debris orbits at the various 
Space Station altitudes is the cause of the scattered effect in the data of 
Figures 10 and 11. The "spikes" in the data of Figure 12 at certain Space 
Station inclinations are caused by the possibility (for some nodal orien- 
tations) for the Space Station to be coplaner with the debris satellite but 
rotating in the opposite direction. This situation manifests itself in equation 
(4) as an unusually large value for d. 
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Figure 12. Total collision probability as a function of Space Station orbit inclination. 




Probobility of Collision 













Probobi lity of Collision 



Figure 14. Probability of collision with individual debris objects 
as a function of apogee altitude of the debris orbit. 
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Figure 15. Variation of collision probability for an individual debris 
object as perigee altitude is varied through the target altitude. 






The effect of the out-of -plane angle, 5 , on the collision probability is 
illustrated in figures 16 and 17 for a few typical debris orbits. These curves 
show that the collision probability is greatest when the debris orbit is coplaner 
with the target orbit. As 5 approaches zero degree a somewhat abrupt increase 
in collision probability occurs compared to a gradual increase as 6 approaches 
180 deg (retrograde coplaner orbits). The major reason for this trend is the 
orientation of the relative velocity vector as a function of 6. In general, high 
probabilities result when the relative velocity vector (and thus the trajectory 
of the debris object) is almost parallel to the direction of travel of the target 
satellite; that is, when its direction cosine is a maximum. Since the inertial 
velocities of the two objects are generally of the same order of magnitude, a 
small change in 5 from zero degree immediately causes the relative velocity 
vector (roughly - V^) to assume a large angle with the path of the target 

vehicle. This acute angle gradually decreases as 6 varies to 180 deg (the 
retrograde position) where the relative velocity vector is again nearly parallel 
to the direction of travel of the target vehicle. 

The collision probability was parameterized over miss distance 
(Fig. 18) and mission duration (Fig. 19). These curves emphasize that 
using large vehicles in earth orbital missions of long duration result in high 
collision probabilities. In Figure 18, the linearity and slope of the logarith- 
mic plot at lower probabilities (less than 0.5) show that the collision probabi- 
ity is approximately proportional to the square of the miss distance. For 
example, if the miss distance is doubled from 50 m to 100 m, the collision 
probability is increased approximately four times. This curve indicates 
that it is almost a certainty that a trackable object will pass within 1/2 km 
of the Space Station at least once during the 10-yr mission. Figure 19 shows 
that the collision probability is linear with mission duration. The circled 
point on Figures 18 and 19 represents the proposed Space Station mission. 


CONCLUSIONS 


If the present trend in satellite population growth continues, the colli- 
sion hazard to men and equipment in earth orbit will become more and more 
severe. It seems inevitable that eventually the development and utilization of 
a collision avoidance system will become mandatory. Indeed, collision 
probabilities for the proposed modular Space Station of 0.01 to 0.04 (for the 
altitudes being considered) may indicate an immediate need for such a system. 
Some alternatives (or additions) to the above remedy are (1) stricter controls 
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Figure 16. Effect of the out-of-plane angle, 6 , on collision probability 
for some typical debris orbits holding apogee altitude constant. 
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Figure 17. Effect of the out-of-plane angle, 5 , on collision probability 
for some typical debris orbits holding perigee altitude constant. 





Probo bi lit y of Collision 



Figure 18. Total collision probability for Space Station as a 
function of miss distance (target sphere radius) . 
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function of mission duration. 



on the amount of debris placed in orbit from a single launch, (2) the reser- 
vation of certain altitudes for manned vehicles or other satellites requiring 
collision protection, and (3) restrictions on the rate of launching objects 
into earth orbits (nationally and internationally). 
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APPENDIX A 


• • # 

DERIVATION OF RELATIVE VELOCITY COMPONENTS X, Y, Z 
AND ORIENTATION ANGLES OF DEBRIS PATH e AND <s 


Suppose the origin of a nonrotating right-handed system of axes is 
placed at the center of the earth and is defined to be an inertial reference 
frame. The origin of the relative coordinate system (Fig. 2) moves at a 
velocity (with respect to the inertial frame) equal to the velocity of the target 
sphere . The immediate objective is to transform the inertial velocity 
vector of the debris satellite into the moving coordinate system. 


The x-axis of the relative coordinate system was specified in the text 
as being a curvilinear axis along the path of the target vehicle. To simplify 
the vector transformation, the curvilinear system can be replaced by a rec- 
tangular Cartesian coordinate system where the y and z axes remain the same 
but the x_axis is tangent to the path of the target satellite. This allows the 
use of the standard transformation between two rectangular Cartesian coor- 
dinate systems, one fixed and one moving. Very little accuracy is sacrificed 
by doing this since computations of the relative velocity components are 
required only in the vicinity where a collision is possible, that is, near the 


origin (within the interval - — ^ x ^ — ) of the relative coordinate 

z z 

system. Under this constraint, the difference between the Cartesian system 
and the curvilinear system is negligible. 


By using the standard transformation between the two coordinate 
systems, the velocity of the debris satellite in the inertial system can be 
written as 

v"=V~ + V~ + "c 17 x p (A- 1 ) 

d r t t 


where 

V = the velocity of the debris object in the relative coordinate system. 

r 

p = the position vector of the debris object in the relative coordinate system. 

u = the rate of rotation of the relative axes with respect to the inertial system. 

The required vector is V . Solving equation (A-l) for this vector yields 

r 


31 



V 

r 


x 


P 


( A-2) 


= V d - V t 


co. 


Writing this equation in matrix form in terms of the three-dimensional com- 
ponents of each term results in 



(A-3) 


Figure A-l shows the orientation of the inertial velocity vector of the debris 
object in terms of the flight path angle, y , and the angle between the debris 
and target orbit planes, 6 (6 is derived in Appendix G) . Resolving this 
vector into its components in the directions of the relative coordinate axes 
yields 


- -V, cos y cos 
d 

= V , sin y 
d 

= V , cos y sin 5 
d 


6 


(A-4) 


Since V remains along the x-axis of the relative coordinate system, 
its three components are 


x 

t 


t 

t 


-V 

t 

0 

0 


(A-5) 


Because of the way the relative coordinate system was defined (p. 4 
second paragraph) it undergoes rotation about the z_axis only. Therefore, 


= 0 (A-6) 

oo =0 

y 

00 -00 
Z t 
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The rotational rate jj is computed as 

2n 

iO 9 

t T 


where 


T = the orbital period of the target satellite. 

The relative velocity components of the debris satellite are computed 
only when the debris satellite passes through the orbit plane of the target 
satellite. Therefore, the last matrix in equation (A-3) becomes 


P z - 0 . 

Substituting (A-4), (A-5), (A-6), and (A-7) into (A-3) yields 


X 


-V _ cosy cos 5 


’ -V 


0 


X 



d r 


t 




d 

y 

= 

V , sin y 
d 

- 

0 

- 

0 

X 

y d 

z 


V . cos y sin 6 
d 


0 


00 

t 


_ 0 


(A-7) 


(A-8) 


Carrying out the matrix operations yields the relative velocity components 

x = V - V, cos y cos 6 + co y , . 

t d r t y d ( A-9) 

y = V sin y - go , x , 
d t d 


z = V cos y sin 6 
d 

The product V cosy in the x and z equations can be replaced by 


C /R . where 
1 d 


V R d v d cos Y ' 

is the angular momentum associated with the orbit of the debris satellite. Also 
y can be written as a function of the true anomaly of the planar intersection 
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line by noting that y is approximately equal to R^. Therefore, 

. _ d , v _ _d_ / p \ ep sin ryg ^ (A-10) 

y_ dt K d^ dt V 1 + e cos 77 / ( 1 + e cost ]) 2 

2 2 

Multiplying the right side of equation (A-10) by /R and noting that 

C 1 = R d^ ’ 
equation (A-10) changes to 

* - pe sin 77 (A-n) 

y ( 1 + e cost)) ,<! R^ 

Since 


R‘ 


d ( 1 + e cos rfp 


equation (A-ll) reduces to 


C«e 

y = — — sin 77 


( A-12) 


Since 


p 2 

,= 

/j.e 


— sin i) 

'i 


(A-13) 


The relative velocity components are finally written as 


x = V - — i 
t R , 


cos 


6 + 




fd 


( A.— 14) 


. jue 

y=— sin 17 


z = — sin 6 
R d 
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(A-15) 


Referring to Figure A-2, the debris path orientation angles are 
0= tan’' (|) 

and 

f tan "' (i) ■ 

Comparing equations (A-15) with equations (A-9), it is found that 
/3=f( x d ) . (A-16) 

The technique used to calculate the collision distance (Appendix B) dictates 

that/3 remain constant as x, varies. Therefore, the ui x. term in equations 

d t d 

(A-9) (the y equation) is considered to be zero since its value will always be 
small (less than 1 m/s) . 
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APPENDIX B 


DERIVATION OF THE COLLISION DISTANCE, d 


As stated in the text, the only varying quantity associated with succes- 
sive passes of the debris satellite through the orbit plane of the target satel- 
lite is the horizontal position of the intersection point. Figure 6 shows 
several of these successive passes. This figure shows that if the debris 
satellite passes through the target orbit plane within a distance of d/2 from 
the y axis, then the debris object will intersect the sphere, thus defining a 
collision. The distance d is a function of the three-dimensional direction of 

travel of the debris object and the altitude difference at plane intersection, y . 

d 

In order to determine d, the equation of the target sphere is first 
solved simultaneously with the equations of the line of travel of the debris 
object. This yields the coordinates of the points at which the debris path 
intersects the sphere. The position of the debris intersection with the x, y 
plane ( x^, y ) is then varied in the x direction to generate the range of values 

of x^ which produces real ( mathematically) intersection points with the target 

sphere. The altitude difference y^ is maintained constant as are the |3 and 9 

angles defining the direction of travel. This range of values of x is equal to 

d 

d. The derivation follows: 

The equation of the sphere which represents the target satellite is 

x 2 + y 2 + z 2 = R 2 , (B-l) 

where R is the miss distance ( or sphere radius) . The equations which specify 
the debris satellite trajectory, assuming that it travels in a straight line and 
has the same direction as the relative velocity vector in the relative coordinate 


system, are 


Y= ” X d^ tan/3+ y d 

(B-2) 

z = (x _ x ,) tan 6 
d 

(B-3) 
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The parameters and are the x and y coordinates, respectively, of the 

point v/here the debris satellite passes through the orbit plane of the target 
satellite. Figure 6 shows this geometrically. To find the intersectio n points 
of the: debris trajectory with the target sphere, equations (B-2) and ( B-3) 
are substituted into equations(B-l) and the terms are collected on x 2 and x, 
yielding 

a t x 2 + bjX + Cj = 0 (B-^ 

wlvere 

a^ = 1 + tan 2 /? + tan 2 0 (B _ 5 ) 

b 1 = 2y^ tan /? - 2x^ ( tan 2 /? + tan 2 0) 

Ch = x 2 , tan 2 /? - 2x ,y. tan /? + y 2 + x 2 tan 2 0 - R 2 
1 d d d J d d 

Solving equation (B-4) for x, the x-coordinates of the points where the li n e 
of travel cuts the sphere are obtained: 

bj. ± n/ b t 2 - 4a t ^ (B_6) 

2aj 2aj 


The y and z coordinates are given by equations (B-2) and (B-3) . 


In equation (B-6) above, if b t 2 = 4a 1 c 1 , the debris trajectory is tan- 
gent to the sphere; if b t 2 > 4sL t c t , it intersects the sphere at two poin' cs; if 
bj 2 < 4ajCi, no intersection occurs. It is required to find the range of values 
of x^ which satisfy the inequality 


bj 2 s 4a 1 c 1 


(B-7) 


Substituting equations (B-5) into equation (B-7) and collecting te rms on x^ 2 
and yields 


ax 
2 d 


b x _ 
2 d 


+ c ^ 0 


(B-8) 
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where 


5i = -4( tan 2 /3 + tan 2 0) 

V 8y d tan fi 

c 2 = 4[(R 2 - y 2 ) sec : 4 + R 2 tan 2 /?] 
1 'he solution of the above inequality is 




The col lision distance, then, is 



( B-9) 


(B-10) 


( B-ll) 


This is th e interval along the x-axis through which the debris object can pass 
and inters ect the sphere. 
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APPENDIX C 


AREA DETERMINATION OF PROBABILITY VERSUS w HUMP. 


To determine the dimensions of the humps in Figure 8, it is noted 

that the maximum probability (for each A£2 increment) as a function of a o 

occurs when the altitudes of the two bodies are equal at the line of intersection 

between the orbit planes. When this occurs the value of y goes to zero and 

d 

the cap that is cut off the sphere in Figure 7, becomes a hemisphere 

causing the collision distance d to reach a maximum. The true anomaly of 

the intersection line ( measured from the line of perigee in the debris orbit) 

at which R = R is calculated as 
d t 

-1 

t] = cos 

max. 


P - R t 

eR. 


(C-l) 


The sum of the true anomaly, r\, and the argument of perigee, cu , is the 
argument of latitude, u. Thus, 


CO = U - T1 

max. max. 


(C-2) 


Substituting the above value of co m into the probability expression 
[equation (5) in the text] yields the maximum collision probability as 

P = d ( W max.’ . (C-3) 

max. — 


where 

l?max = the height of the probability hump of Figure 8. 

To obtain the width of the hump, A co , it is required to find the values 
of co where the probability goes to zero. The probability expression is set 
equal to zero yielding 

P, = d/AS = 0 
lap 
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or 


d = 0 

From Appendix B, 



®2 


Therefore to satisfy equation (C-5) 
b 2 ^ - 4a 2 c 2 =0 , 

where 

a 2 = -4 ( tan 2 /3 + tan 2 0) 
b 2 = 8y d tan/3 

Cg= 4[(R 2 _ y^ 2 ) sec 2 0 + R 2 tan 2 /3] 


(C-5) 

(C-6) 

(C-7) 

(C-8) 


Substituting equations (A-15) and (C-8) into (C-7) and solving for the 
altitude difference, y , yields 



where y and z are derived in Appendix A. The quantities y and y are the 

dj d 2 

altitude differences between the two bodies where the collision probability 
goes to zero. These two values are then used to calculate the orbital radius 
of the debris satellite ( at the line of intersection between the orbit planes) 
at which the collision probability goes to zero as 





(C-10) 


The position ( true anomaly) in the orbit of the debris satellite where this 
radius occurs is 



(C-U) 
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where 


p= the semilatus rectum 
e = the eccentricity. 

The arguments of perigee corresponding tor]* and rj 2 are expressed as 

w i,2= u “^1,2 » ( C-12) 

where 

u = the argument of latitude (in the debris orbit) of the intersection 

line. 

The difference between uo { and co 2 is the width of the probability versus 
hump. Thus, 

Aw = w 2 - W! = (u - t] 2 ) - (u -rfj) 


or 


Aw = Tjj -T ) 2 • (C-13) 

Aw is typically less than . 5 deg for a miss distance of 50 m. The area of 
the probability versus w hump is then approximated as one-half of an ellipse 
with the minor axis equal to Aw and the maior axis equal to 2 P 

max. 

(Appendix D for justification) . The approximate area of one hump is 

A = J tt P Aw . (C-14) 

max. 
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APPENDIX D 


MATHEMATICAL JUSTIFICATION FOR 
ELLIPTICAL APPROXIMATION OF INTEGRAL 


Figure D_1 shows the target sphere and the two positions at which 
the debris path can be tangent to the sphere. These two tangent points lie 
on a small circle of the sphere. An edgewise view of this circle is shown 
in Figure D-2 as seen by an observer looking down the x_axis in the +x direc- 
tion. The plane of this circle is always perpendicular to the y-z plane since 
y^, 6 , and / 3 remain constant for successive passes of the debris satellite 

through the target orbit plane and the only variable is x^. In figure D_2, the 
radius of the small circle is r and 

r 2 = R 2 _ (y cos o' ) 2 . (D-l) 

Since 

« = tan ~ 1 ('z' 



then 


cos a 


z 

\T y 2 + z 2 


Equation (D-l) can then be written 



( 


if \ 

y z + z 2 / 


(D_2) 


From Figure D-l, d can be related to r in terms of the relative velo- 
city components by taking the dot product of and f, a unit vector in the 

x-direction, as follows: 

V • f= V cos e = x . (D-3) 

r r ' 


Therefore 


cos e = 


x 


V 

r 


(D-4) 
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Target Sphere 


To Earth 

Figure D-2. Edgewise view of the small circle in the target sphere 
defined by passage of the debris object through the target sphere. 




Squaring equation (D-4) and substituting it into the trigonometric identity 


sin 2 e = 1 - cos 2 e 


(D_5) 


yields 


. ? x- 

sin € = 1 - — 2 

r 


y 2 + z 2 


v 


• 9 • ? 

+ z u 

- -rrj m ry 

x + y + z 


(D_6) 


From Figure D-l 
2r 

sin e = — r 


(D-7) 


Squaring equation (D-7), substituting equation (D-6), and solving for d 2 
results in 


d 2 = 4r 2 


x 2 + y 2 + z 2 

ry ry 

y z + z 


(D-8) 


Substituting (D-2) into (D-8) gives 

_> 2 
■ y‘ 


d 2 = 4 


r2 - y d 2 (ftf) 


x 2 + y 2 + z 2 
y 2 + z L ~~ 


(D-9) 


Carrying out the multiplication yields 


■2 _ 4R 2 ( x 2 + y 2 + z 2 ) _ 4z 2 (x 2 + y 2 + z 2 ) , 

( y 2 + W 


y 2 + z 2 


and upon rearranging, 


d 2 + 


2 ( '2 


4z 2 (x 


v + z 


z 2 ) 2 _ 4R 2 (x 2 + y 2 + z 2 ) 


( y 2 + z 2 ) 2 


y T TT r 


(D-10) 
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Both sides of equation (D-10) are divided by the right side to give 

X, 2 


d 2 

« I x 2 + y 2 + z 2 

4R 2 1 iTiV - 


rx i+ y‘+z‘i t> 2 r y 2 + z' z -i 

L F~rz J J 


= 1 


(D-ll) 


Since 


lap AS 


then 


d 2 = P, 2 AS 2 
lap 


(D-12) 


Substituting (D-12) into (D-ll) yields 

r> 2 
lap 


4R 

AS 


2 r x 2 + y 2 + z 2 q 

i» L J 




= 1 


(D-13) 


The altitude difference, y , can be expressed as a linear function of 

d 

rj (the true anomaly at which the intersection line occurs) by writing the 
expression for y as a Taylor series. The expression for y^ is 


y^_ 


d = f (tj) = - 


R 


+ e cos 7) 


(D-14) 


The Taylor series of f(T)) expanded about the point tj is 

f (rj ) = f (tj )+ f'(rj ) (tj - r? ) + f (v -V ) 2 + ••• 

o o 2! 


(D-15) 


Substituting (D-14) into (D-15) and neglecting higher order terms (for small 
At) ' s) yields 

P( 1 - e sin t) ) 


y d = 


1 + e cos r] 


- R 1 - [ P(l -esfan o ) 1 

1 J L ( 1 + e cos t) ) 2 J 
o 


At) 


(D-16) 


where 


At) = tj - T] 
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This linear function is a close approximation to the actual function if A?] is 
sufficiently small. A typical value for Atj is less than .5 deg. 


Letting the constants in equation (D-16) be 
P 


a = 


1 + e cos r] 


- R. 


and 


b = 


-P( 1 - e sin rj ) 
( 1 + e cost] ) 2 


(D-17) 


(D-18) 


equation (D-16) then becomes 

y = a + b Arj . (D-19) 

If is chosen such that a= 0 [equation (D_17)J , then 

y d = b A?7 . (D-20) 

Substituting equation (D-20) into equation (D-13) yields 


P 2 

lap An 2 

' \ T ' ■, ^ I r? r ir“ _L 


4R 2 

• 7 * [T I 'J 

" X + jr + z “ 

+ R 5 * p jr + z z " 

= 1 

• 

(D-21) 

AS 2 

• y 

Y + 7T 

|_ z z 





This equation shows that the probability of collision versus rj curve is ellipti- 
cally shaped under the assumptions that (a) the velocity components of the 
debris satellite remain constant over the width (At]) of the ellipse, and (b) 
the interval At], for which a non-zero probability exists, is sufficiently small 
to insure that y^ is a linear function of t] . These assumptions are not true only 

when the debris satellite radius of perigee is within several target sphere 
radii, R, of the radius of the target satellites orbital radius. 
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APPENDIX E 


EXTENSION OF PROBABILITY PER LAP EXPRESSION 
TO PROBAB ILITY FOR THE ENTIRE MISSION 


Let P be the probability, for each lap, that a debris satellite will 

i£Lp 

collide with the target satellite. The probability, then, that a collision will 

not occur during the first lap of the mission is 1-P, . The theorem of 

lap 

compound probability states that if A and B are any two independent events, 
then 


p(AB) = P ( A) p(B) . ( E-l) 

Stating equation (E-l) informally, the probability that both A and B happen is 
the probability that A happens times the probability that B then happens. Now 
associating A with the probability of noncollision for the first lap and asso- 
ciating B with the probability of noncollision for the second lap, the probability 
of noncollision for both laps is 


P 

nc 


(1 - 




(1 


P lap 2 ^ 


or 


p = ( i _ p ) 2 

nc lap 


(E-2) 


since the probabilities of noncollision are equal for both laps (P = P ) 
The probability of noncollision for L laps is then ^ 

P = (1 - P. ) L . (E-3) 

nc ' lap' v ’ 


Since 


^collision + P noncollision ^ 


( E-4) 


then the probability of collision with the jth satellite for a mission duration 
consisting of L laps is 


P 

m. 

J 


(1 - 


P lap) L 

J 


(E-5) 
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APPENDIX F 


COLLISION PROBABILITY BETWEEN THE TARGET SATELLITE 
AND ANY ONE OF THE TOTAL NUMBER OF DEBRIS SATELLITES 


Let be the probability of collision between the target satellite and 

j 

the jth debris satellite for the total mission. The objective is to derive an 
expression for the probability that a collision will occur with any debris object 
considering the entire debris population. The principle of compound probability 
as explained in Appendix E applies in this case also. The difference here is 
that the individual satellite probabilities, P , will not be equal in general. 

j 

Thus, the total collision probability is written as 


d 

P TOT 1 " JT ( 1 ” P m ^ 

J=1 j 

where 


n^ = the total number of debris satellites in earth orbit. 
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APPENDIX G 


DETERMINATION OF THE ANGLE BETWEEN TWO ORBIT PLANES 


Figure G-l shows the projections of two orbit paths onto a spherical 
earth. The orbits represented by these projections are inclined to the equa- 
torial plane by amounts i and i . The angle between the planes, 6, can be 

d t 

found in terms of the known quantities, i^, i^, and A ft by considering the 

oblique spherical triangle ABC. The relations between the sides and angles of 
this triangle are completely specified by NapieFs analogies. The equations 
are as follows: 


sm id; - t d ) 


sin 2 (ij. 

+ *d) 

cos 2 ( i^ 

- i d ) 

cos 2 ( i' 

+ i d ) 

sin \ ( <j 

- A) 

sin 2 (a 

+ a) 

cos \ ( cr 

— i ; 

- A) 
\ 


COS 


tan 

2 ( V - 

A) 

tan 

2 ( A £2 ) 


tan 

h(cr + 

A) 

tan 

2 ( A £2 ) 


tan 

id; - 

i d ) 


cot 4 6 


tan 

id; ♦ 
1 

i d ) 


cos ( 6 


(G-l) 


( G-2) 


(G-3) 


(G-4) 
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Projection of Target Orbit 
Onto the Earth 
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Orbit onto The Earth 

Figure G-l. Geometrical description of the mutual orbit orientation 
between the debris and target orbits. 


Equations (G-l), (G-2), and (G-4) are sufficient to find <5 

(o "— ia \ 

solve equations (G-l) and (G-2) for 

obtain 


First, 

(ct-\) (tr+X) .. , , 

- and — - — respec tively, to 


= tan' 

2 


sin \ (i t - i d ) esc h (i t + i d ) tan \ AfiJ (G-5) 


and 


— * — = tan 1 |^cosi (i^ - i d ) sec \ (ij + i^) 

tan \ A £2 J . (G- 6 ) 

Substituting equations (G-5) and (G- 6 ) into (G-4) to elimin ate the unknown 
quantities and then solving for 6 yields 

6=2 cot -1 £ tan \ (i* + i d ) cos tan -1 (j) sec ta.n _1 (k)J , (G-7) 

where 

j = cos \ ( i’ - i d ) sec \ ( i’ + i d ) tan \ A S 2 


and 

k = sin \ (i’ - i d ) cse |(i^ + i d ) tan \ Afi 

Noting from Figure G-l that i’ = ir - i^, and substituting this into equation 
(G-7) results in 

6=2 cot -1 ^cot 5 (i^ - i d ) cos tan 1 (j T ) sec tan -1 (k f ) J , (G- 8 ) 

where 

j’ = sin \ (i^ + i d ) esc \ (i^ - i d ) tan \ A £2 

and 

k’ = cos I (i t + i d ) sec i (i t - i^) t;rn \ A £2 
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